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Abstract 

A theoretical investigation is made of the dispersion characteristics of plasmons in a two- 
dimensional periodic system of semiconductor (dielectric) cylinders embedded in a dielectric 
(semiconductor) background. We consider both square and hexagonal arrangements and 
calculate extensive band structures for plasmons using a plane-wave method within the 
framework of a local theory. It is found that such a system of semiconductor-dielectric 
composite can give rise to huge full band gaps (with a gap to midgap ratio ~ 2) within 
which plasmon propagation is forbidden. The most interesting aspect of this investigation 
is the huge lowest gap occurring below a threshold frequency and extending up to zero. The 
maximum magnitude of this gap is defined by the plasmon frequency of the inclusions or the 
background as the case may be. In general we find that greater the dielectric (and plasmon 
frequency) mismatch, the larger this lowest band-gap. Whether or not some higher energy 
gaps appear, the lowest gap is always seen to exist over the whole range of filling fraction 
in both geometries. Just like photonic and phononic band-gap crystals, semiconducting 
band-gap crystals should have important consequences for designing useful semiconductor 
devices in solid state plasmas. 



I. INTRODUCTION 

Ever since the discovery of quantum Hall effect, the investigation of electronic transport in the 
low-dimensional systems has witnessed many important and fascinating results, both theoretical 
and experimental^. The ever increasing interest in these systems is attributed partly to advance- 
ment in the modern semiconductor technology (e.g., molecular beam eptaxy and electronic beam 
lithography) that makes it possible to fabricate the nanostructures with controllable chemical 
composition and geometric structures^ — quantum-wells, -wires, -dots, -rings, -crossbars, -balls 
are a few to name. The keystone behind the intense research interest in such systems world-wide 
is esentially the fact that their dimension(s) are (or have to be) comparable to the de Broglie 
wavelength, so as to observe the quantum-size effects. 

Many novel discoveries in physics have originated from the study of wave propagation in 
periodic structures. Examples include x-ray and electron diffraction by crystals, electronic band 
structure, Bloch oscillator, and holography. Analogies between subfields of physics have also 
opened amazingly new fruitful avenues in research. An exciting example is the recent discovery^ 
of periodic dielectric structures that exhibit what is called "photonic band gap", by analogy 
with the electronic band gap in semiconductor crystals. Photons within the frequency range of 
the gaps are completely forbidden, so that atoms in such materials are denied to spontaneously 
absorb and re-emit light in this region; this has obvious implications for highly efficient lasers. 
Since the interesting phenomena emerging from photonic band gap structures, such as the in- 
hibition/prohibition of spontaneous emission, the formation of localized states of light, and the 
photon-atom bound states, are all the consequences of the existence of the gap(s), a major part 
of research efforts has concentrated on the search for such photonic crystals. Within almost a 
decade the field has prospered much more than Eli Yablonovitch and Sajeev John could have 
imagined^. 

It did not take long before photonic crystals aroused interest in the analogous studies in other 
systems as well; we refer to the ones allowing, for example, elastic/acoustic waves^, and spin 
waves^. The composite sytems exhibiting complete or pseudogaps, within which the respective 
waves are forbidden, are useful both from the practical and fundamental point of view. From 
the practical point of view, such systems can be used to design the filters that prohibit the 
propagation of respective waves at certain frequencies while allowing free propagation at others. 



From the fundamental point of view, the propagation of such waves in a shghtly disordered 
composite system can lead to the novel phenomenon concerned with the Anderson localization. It 
is noteworthy that the authentic seeds of classical wave localization were disseminated in a classic 
paper by Sajeev John in 1984^. A systematic account of the theoretical and experimental work 
on electronic-, phononic-, photonic-, and vibrational-band gap crystals is given in an extensive 
review article by the author^. 

In this work we are concerned with a two-dimensional (2D) periodic system of semiconducting 
(dielectric) cylinders embedded in a dielectric (semiconductor) background. In fact, it could be 
considered a bulk-mass of semiconductor with periodically punched infinitely long thin holes. 
The cross-sections of these cylinders (or punched holes) are considered to be much less than the 
de Broglie wavelength, so as to neglect the quantum-size effects. We believe that such geometric 
structures should be easily realizable through the aforesaid nanostructure growth techniques. 
Our principal interest is the prospect of achieving a complete plasmonic band gap (PBG); this is 
defined to be a stop band within which the plasmon propagation is prohibited irrespective of the 
direction of propagation and of the polarization of the plasma wave. The existence of PBG in 
such composites should have important consequences for designing filters for the plasma waves 
and for the fundamental physics interest associated with the Anderson localization thereof. In 
fact, the possibility of experimental realization of the Anderson localization of plasma waves in 
random and incommensurate layered systems was predicted almost a decade ago^. 

It is noteworthy that our model theory is based on the full vectorial Maxwell's equations. We 
also consider an applied magnetostatic filed oriented parallel to the axis of the cylinders (which 
is taken to be the z-axis). In order to calculate the band structure of such a binary composite, 
we employ a plane-wave method and Fourier-series-expansion of the position-dependent material 
parameters. It should thus be pointed out that we do not make use of any messy electromagnetic 
boundary conditions usually employed for inhomogeneous systems. The semiconductor medium 
in the inclusions and/or background is characterized by the frequency-dependent (local) dielectric 
function. The first hand effort of using local theory for photonic crystals has been made in the 
recent past^. We neglect the damping effects and hence ignore the absorption. In the state-of- 
the-art high quality semiconductor systems available, this is thought to be quite a reasonable 
approximation. Thus our purpose is to calculate the band structure of a neat and clean 2D 
periodic binary semiconductor-dielectric composites. 



A number of investigators have calculated such dispersion characteristics of ID periodic sys- 
tems (superlattices, for example) both with and without an applied magnetic field^°^^^. However, 
no one, to our knowledge, has ever embarked on the possibility of observing complete band gaps in 
the systems considered. The authors were, most of the time, interested in reporting the physical 
conditions leading to the existence of the bulk bands and the surface plasmons or magnetoplas- 
mons in the infinite and semi-infinite superlattice systems. This seems to be true irrespective 
of the systems (compositional or doping (n-i-p-i) superlattices) considered, methodology (classi- 
cal or quantal) employed, and the orientation of the magnetic field (Voigt-, perpendicular-, or 
Faraday-geometry) pondered. The same remark is valid for the work on the ID lateral (quantum 
wire) superlattices"*^^"^^. 

The rest of the paper is designed as follows. In Sec. II we derive the secular equations 
required to compute the plasmon band structures of the 2D binary semiconductor composites. 
The extensive numerical computational results are discussed in Sec. III. This includes results 
both on square and hexagonal arrangements. A brief summary in Sec. IV concludes our finding 
in this paper. 

II. METHODOLOGICAL DETAILS 

We consider a 2D periodic system of semiconducting cylinders embedded in a different semi- 
conducting background. For the sake of generality, we consider both cylindrical inclusions and 
the background materials to be semiconductors to be characterized by the frequency-dependent 
(local) dielectric functions. An external magnetostatic field {Bo) is assumed to be oriented 
along the axis (z) of the cylinders. We are interested in the nonmagnetic materials, which im- 
plies that B = H in the Maxwell equations. In the situation at hand {Bo \\ -2),the dielectric 
tensor e is simplified by the symmetry requirements such that: exx = ^yy, ^yx = —^xy, and 

^zx ^xz ^zy ^yz U. 

We start with full vectorial Maxwell's curl field equations. After eliminating the magnetic 
field variable B, we obtain the following wave field equation for the macroscopic electric field 
vector (i?): 

V x{V X E)-qli-E = (1) 



where e is a 3 x 3 matrix and go = ^/c is the vacuum vector; c being the speed of hght in vacuum 
and 00 is the angular frequency of the wave. Eq.(l) can be written in the matrix form as follows: 
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Here dz = d/dz-, etc. Eq. (2) can be solved for three different situations: (i) kz = 0, and k^\ ^ 0; 
(ii) kz ^ and k^\ = 0; and (iii) kz ^ and /cy ^ 0. We are here interested in the first case 
— with kz = but /cji 7^ 0. Here k\\ is the 2D Bloch vector and k^ is the z-component of the 
propagation vector. 

As such, eq. (2) splits in two independent sets of equations (with kz = 0): 



V^Ez + qlezzEz = 



(3) 



and 



{ql^xx + dl)Ex + {qlexy - dxdy)Ey = 
-ql^xy - dxdy)Ex + {qle^x + dl)Ey = 



(4.1) 
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where eij are all functions of position. Evidently, eq. (3) represents the propagation of the plasma 
waves which remain independent of the magnetic field; whereas the set of eqs. (4) refers to the 
magnetoplasma waves. We will therefore derive two independent secular equations corresponding 
to the two situations. It is, however, noteworthy that both equations (3) and (4) are valid for 
arbitrary inhomogeneity. 



2.1. Magnetic-filed-independent waves 

Here we are concerned with the plasma waves governed by eq. (3). We focus our attention on 
an inhomogeneous medium which exhibits spatial periodicity. That means that the position- 
dependent material parameters tzzi'f^) may be expanded in the Fourier series. We write 



(^zz{r) = J2^zziG)e 

G 



iG-r 



(5) 



where G is the reciprocal lattice vector that has the corresponding dimensionality. The z- 
component of the field vector must satisfy the Bloch theorem: 



i?,(f,t)=e*(^-^^-'^*)^i?,(G)e^ 
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(6) 



Equations (5) and (6) must be substituted in eq. (3). A straightforward algebra, then multiplying 
the whole equation by exp{—iG -r), and integrating over the unit cell yields the following result: 
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For binary composites every unit cell is composed of only two materials — inclusions (i) and 
background (b) characterized by e!,^ and e^^. In addition, we must specify the filling fraction / 
and (1 — /), respectively, for materials i and b. We consider ezz to be given by 
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where Ac is the area of the unit cell and the integration is over the unit cell. A simple conceptual 
aspect of the system leads us to define 
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where the structure factor F{G) is defined as 
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where Ji{x) is the Bessel function of the first kind of order one, the integration is limited to a 
cylinder of radius Tq, and the filling fraction 
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nij-Q/aY , for square lattice 



:ii: 



27r 



/^(ro/a)^ , for hexagonal lattice 

where a is the period (or lattice constant) of the system. With the aid of eq. (9), one can cast 
eq. (7) in the form (see Appendix): 
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where 
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where e^^ and ujpj are, respectively, the background dielectric constant and the screened plasma 
frequency of the medium concerned; j = i, b. 

Equation (12) can be cast in the form of a standard eigenvalue problem: 

9' 

where the dimensionless Q = uj{a/2TTc), k = K[a/27r), and g = G{a/27i). The matrix A is defined 
as follows: 

A=B C (15) 

and the matrix elements of B and C are given by 

Bm' = Ph,9' + RF{g-^){l - %gO (16.1) 
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where 



Q' = Q{a/2ncy = e^^U + e^fij,(l - /) (17.1) 

S' = S{a/27rcf = e^^l, - e^^l, (17.2) 

where Vtpj = ujpj{a/2TTc) is the normalized screened plasma frequency. Equation (14) is, in fact, 
performed at the computational level for both the geometries considered in the present work. 
The existence of matrix A implies that its elements must be computed for every value of k, rather 
than for every value of k and Q as implied by eq. (12). A drastic saving in computational time 
is achieved by transforming eq. (12) into eq. (14). The penality one has to pay for this benefit 
is that one is required to invert matrix B and perform the multiplication B C- 



2.2. Magnetic-field-dependent waves 



Now we turn to the magnetoplasma waves governed by eqs. (4). Invoking the Bloch theorem for 
the field components {E^ and Ey), making Fourier expansion of the position-dependent material 
parameters {e^^ and e^y), and normalising the quantities, just as before, leads us to write eqs. 
(4) in the form 
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where functional A{g) and B{g) are defined as 
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where k and g are just as defined before and 
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with sub/superscript j = i, h. Here Vt^j = {oOcja/l-Kc) is the normalized electron cyclotron 
frequency (see Appendix). 

Let us now rewrite eqs. (18) symbolically as follows: 



A,E,{^)+A^Ey{g')=Q 
A^E^{^)+A^Ey{g')={) 



(22.1) 
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where the matrices Aj;j = 1,2,3,4; are defined such that the matrix elements thereof are given 

by 

Ai(g,g') = A{g - g) - {k + g)y{k ~ g)ySg'^g (23.1) 

A2iM') = B{g-g') + {k + g').{k - g')y6g,^g (23.2) 

^3(5,,-) = -B{g-g') + {k + g').{k-g')y6g,.g (23.3) 

^4{s,a') = A{g- g) - {k + g)^ik - g)J^>^^ (23.4) 

Ehminating Ey{g') from eqs. (22) leaves us with 

%''a^ - i^'^s) E,{g') = (24) 

Nontrivial solutions of this equation are given by 



A2 Ai - Ai A. 



(25) 



What is noteworthy is the fact that in this magnetic-field-dependent part of the problem, the 
secular equation cannot be cast in the form of a standard eigenvalue problem. That is because 
the eigenvalues are involved in both the diagonal and nondiagonal elements of the resultant 
matrix. So the only way down to the determination of the eigenvalues is to treat eq. (25) as 
a transcendental function and find its zeros for a given k\\, by inserting an additional loop over 
Q. Naturally, the time consumption during the computation of reliable band structures will be 
enormous. 

It is interesting to note that one can, without much efforts, prove analytically that, in the 
limit of vanishing applied magnetic field (=^ functional B{g) -^ 0), the secular equation (25) 
reduces to 



A-n' I 



(26) 



where / is a unit matrix. Eq. (26) is exactly identical to the secular equation obtained from eq. 

(14). This emboldens our confidence in the adequacy of the formalism presented in this work. 

The last step in presenting the numerical examples is to specify the symmetry of the lattices. 

The specific geometrical arrangements considered in this work are the square and hexagonal 

lattices. Their reciprocal lattice vectors are given by 

2n 
G = — {n^x + nyj)) (27) 

(Jj 



for the square lattice and 

G=— \n^x + ((-n, + 2ny)/^'i)y] (28) 

for the hexagonal lattice. In this work the integers n^ and Uy were limited in the range of 
—10 < {rixjUy) < +10; which implies to 441 plane waves considered. This resulted in a reliably 
very good convergence; at least up to the first fiftieth bands. In this paper we will confine 
ourselves to the waves which turn out to be independent of the applied magnetic field and hence 
are governed by eq. (14). The analytical and numerical details of the magnetic-field dependent 
waves are deferred to a forthcoming publication. 

III. NUMERICAL EXAMPLES 

For the purpose of numerical computation, we specify the semiconductor- dielectric composite to 
be made up of doped GaAs-intrinsic GaAs or doped GaAs-vacuum. This implies that ei = 12.8 
(1.0) for doped or undoped GaAs (vacuum). As stated above, we confinr ourselves to the waves 
which are independent of the applied magnetic field. As such the only other parameter needed to 
accomplish the computation is the dimensionless plasma frequency Qp] two specific values used 
in this work are Qp= 0.5 and 0.75. We will discuss a variety of illustrative numerical results. 
This includes exchanging materials in the cylindrical inclusions and the background medium; in 
addition to the two values of Qp. The whole range of filling fraction is explored in order to discuss 
the regime of at least first ten bands of the respective composites in both square and hexagonal 
lattices (see Fig. 1). 

3.1. Square lattice 

Figure 2 illustrates the band structure (BS) and density of states (DOS) for a square lattice of 
doped GaAs cylinders in vacuum for a filling fraction / = 10%. The plots are rendered in terms 
of the dimensionless frequency Q versus Bloch vector k. The left part of the triptych represents 
the BS in three principal symmetry directions, letting k scan only the periphery of the irreducible 
part of the first Brillouin zone (see the inset on the right). The middle part is the result of an 
extensive scanning of | A; | in the irreducible triangle FXM of the Brillouin zone — the interior of 
this zone and its surface , as well as the principal directions shown in the left part of this figure. 
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Each curve here corresponds to some direction of k. The DOS in the right part of the triptych 
has been calculated on the basis of the scanning in the middle part, which corresponds to 1326 
k points. The three parts of the triptych together demonstrate that there is, indeed, a genuine 
full gap in the frequency range defined by < ^2 < 0.320 and we consider such calculations as 
essential. 

Extensive computation of the BS for the same composite as considered in Fig. 2 reveals that 
one can also achieve some higher energy gaps within the first ten bands for certain values of the 
filling fraction. For instance, for / = 30%, there are three full gaps within the first ten bands — 
the first one is analogous to the one shown in Fig. 2, the second one exists between the third and 
fourth bands, and the third one opens between the ninth and tenth bands. For / = 40%, the 
third gap exists between the sixth and seventh bands, whereas the second gap, although differs in 
magnitude and frequency range, still persists between the third and fourth bands. For / = 50%, 
there are only two gaps — the first one exists below the first band and the second one opens 
between the sixth and seventh. We find that, whether or not the higher energy gaps exist, the 
lowest gap below a threshold frequency Qt (which varies with /) defined as the minimum of the 
first band at the F point, always persists. It has been noted that, for / > 60%, only the lowest 
gap exists. This is true even for the close-packing value (/ = 0.7853) where the cylinders just 
touch each other. 

The dependence of the gap-widths of the lowest (which is also the widest) gap on the filling 
fraction for diversely designed composites in square lattice is summarized in Fig. 3. The curves 
designated as ASj, SAj, ISj, and SIj refer , respectively, to the composites made up of cylindrical 
holes in doped GaAs background, cylindrical doped GaAs in vacuum, cylindrical intrinsic GaAs 
in doped GaAs background, and cylindrical doped GaAs in intrinsic GaAs background; j = 1(2) 
correspond to Qp = 0.5(0.75). The important points one can notice at glance are the following. 
The composites made up of doped semiconductor inclusions in the dielectric background give 
rise to the lowest gap (below Qt) whose width increases with increasing filling fraction (see, for 
example, the curves SAj and SIj; j = 1,2). On the other hand, the composites made up of 
dielectric cylinders in the doped semiconductor background give rise to the lowest gap which 
decreases as the filling fraction increases (see, for example, the curves ASj and ISj). Comparing 
the curves ASj with ISj (or SAj with SIj) leads us to infer that the larger the dielectric mismatch 
(between the inclusions and the background), the wider is the gap-width. Similarly, comparing 
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the respective curves for j=l with those for j=2 reveals that the gap- width becomes larger as the 
plasma frequency increases. It is also noteworthy that the maximum gap-width (of the lowest 
gap) one can achieve is, empirically, defined by < ^2^ < r2p. This implies that the largest gap 
to midgap ratio (GMR) one can achieve is exactly two. This abides by the fact that in order to 
obtain GMR > 1, 1^2 should be > 3i7i; ^2 (^i) being the frequency of the top (bottom) of the 
existing gap. 

3.2. Hexagonal lattice 

Figure 4 ilustrates the band structure and DOS for a hexagonal lattice of doped GaAs cylinders 
in vacuum for a filling fraction / = 10%. The plots are rendered in terms of the dimensionless 
frerquency Q versus Bloch vector k. The left part of the triptych represents the BS in three 
principal symmetry directions, letting k scan only the periphery of the irreducible part of the 
first Brillouin zone (see the inset on the right). The middle part is the result of an extensive 
scanning of | /c | in the irreducible triangle TJX of the Brillouin zone — the interior of this zone 
and its surface , as well as the principal directions shown in the left part of this figure. The DOS 
in the right part of the triptych has been calculated on the basis of the scanning in the middle 
part, which corresponds to 1275 k points. The three parts of the triptych together demonstrate 
that there is a genuine full gap in the frequency range defined by < f] < 0.329. In addition, 
there is another band-gap opening up between the first and second bands — defined by the 
maximum (minimum) of the first (second) band at J (X) point. 

Extensive computation of the BS for the same composite as represented by Fig. 4 reveals that 
there are some other higher energy gaps opening up within the first ten bands for certain values 
of /. For instance, for / = 20%, there are three more full gaps, apart from the lowest one that 
opens below Qt — the first one exists between the first two bands, the second one opens between 
the third and fourth bands, and the third one exists between the eighth and ninth bands. For 
30% < f < 60%, we obtain two more full gaps, apart from the lowest one, with varying frequency 
range and magnitudes. At / = 70%, we left with only two full gaps — the first one is the same 
that exists below the minimum (at F point) of the lowest band and the second one opens between 
the third and fourth bands. For / > 70%, there are no more full gaps except the lowest one. 
The semiconductor-dielectric composites in the hexagonal lattice, just as in the square lattice. 
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are always seen to give rise to the lowest gap persisting below a threshold frequency Qf 

The dependence of the lowest gap-width on the filling fraction for diversely designed com- 
posites in the hexagonal lattice is depicted in Fig. 5. Comparing Fig. 5 with Fig. 3 clearly 
reveals that the gap-widths (of the lowest gap) for the respective composites in the hexagonal 
lattice, over almost the whole common range of filling fractions, are larger than those in the 
square lattice. The rest of the discussion related to Fig. 3 is still seen to be valid. It should be 
pointed out that this trend of achieving larger gaps in the hexagonal lattice as compared to the 
square lattice is the same as has been noted earlier in the cases of photonic as well as phononic 
crystals^. This seems to be physically quite a reasonable result, because the constant energy 
surfaces in the hexagonal lattice are closer to the circular shape than those of a square lattice. 

IV. CONCLUDING REMARKS 

In summary, we have demonstrated that the simple 2D periodic binary semiconductor- dielectric 
composites can give rise to full plasmonic band-gaps. Although we consider an applied magnetic 
field oriented along the axis of the cylindrical inclusions, two independent polarizations of the 
waves have been obtained — Ez \\ z axis and Hz \\ z axis, just as in the absence of an applied 
magnetic field. Only the waves in the latter polarization are magnetic- field dependent. The 
numerical results presented in this work correspond to the former polarization. As such the term 
"full" should be reserved with respect to this polarization. In principle, the band structures 
computed separately for the two polarizations should be superimposed — one over the other — 
to check and claim whether or not the band-gaps are independent of the polarization. This is 
the task aimed at in a future publication. 

Our theoretical results on achieving full band-gaps must be accompanied by a word of warning. 
Apart from the fact that both polarizations need to be considered, we have here also assumed 
the z-component of the wave vector (kz) to be zero. Considering kz finite would give a better 
insight into the problem. We have also ignored the damping effects and the interaction of the 
charge carriers with the optical phonons. The latter could be taken into account by considering 
the frequency-dependence of the background dielectric constants. Additional work addressing 
the effect stemming from the cylinders of different cross-sections and orientation of the external 
magnetic field would be worthwhile. Considering the cylindrical inclusions of finite lengths should 
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be of more practical interest. 

Attractive possibilities exist for the experimental observation of such plasmonic band-gaps, 
both with and without an applied magnetic field, in the semiconductor- dielectric composites as 
investigated here. Resonant Raman scattering, as employed to study plasmons in quantum wire 
superlattices, may be one option worth attempting. 
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APPENDIX 

The dielectric tensor components employed in this work are defined as follows: 
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(29.1) 
(29.2) 
(29.3) 



where e^ is the background dielectric constant, v is the free carrier collision frequency, uOp is the 
screened plasma frequency, and ujc is the electron cyclotron frequency. 

If we also consider the effect of phonons, which, in a way, incorporates the coupling of the 
plasmons (or magnetoplasmons) to the optical phonons, then the background dielectric constant 
e^ has to be replaced by its frequency dependent expression: 



eLl^j 



ujIq — uj — iVuj 



(30) 



where eoo is the high-frequency dielectric constant, F is the optical-phonon damping frequency, 
and Ulo and ujto are, respectively, the longitudinal and transverse optical phonon frequencies 
at the zone center of the Brillouin zone. 
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Figure captions 

Fig. 1. Schematics of a few units cells of square (a) and hexagonal (b) lattices investi- 
gated in this work. The cylindrical inclusions are oriented along the z axis. The 
propagation is confined to the x — y plane. An external magnetic field field is 
orinted parallel to the cylinders. 

Fig. 2. Band structure (BS) and density of states (DOS) of a square lattice of cylindrical 
(GaAs) semiconductor inclusions of circular cross-section in a vacuum background. 
The filling fraction is / = 10%. The triptych is comprised of three parts. In the 
left panel, we plot the BS in three principal symmetry directions, letting k scan 
only the periphery of the irreducible part of the first Brillouin zine (see inset on 
right). The middle panel demonstrates a novel way of plotting the eigenvalues as a 
function of | /c |; i.e., the distance of a point in the irreducible part of the Brillouin 
zone from the F point. The right panel illustrates the DOS. We call attention to 
the full gap below a threshold frequency ^2^=0. 320 and extending down to zero. 

Fig. 3. Dependence of the gap-widths (of the lowest gap) on the filling fraction for di- 
versely designed composites in a square lattice. The curves designated as ASj, SAj, 
ISj, and SIj refer, respectively, to the composites made up of cylindrical holes in 
doped GaAs background, cylindrical doped GaAs in vacuum, cylindrical intrin- 
sic GaAs in doped GaAs background, and cylindrical doped GaAs in intrinsic 
GaAs background. The horizontal dashed lines stand for the plasma frequencies 
VLpi = 0.50 and Qp2 = 0.75. 

Fig. 4. The same as in Fig. 2, but for the hexagonal lattice. The threshold frequency 
for the lowest gap that extends down to zero is Qt = 0.329. Notice a small second 
full gap (hatched region) existing between the first two bands. 

Fig. 5. The same as in Fig. 3, but for the hexagonal lattice. 
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